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1 Introduction 

The purpose of this paper is to discuss analyticity properties of two-point functions 
in the Nambu-Jona-Lasinio model. We point out that the finite cut-off which regulates 
the quark loop can lead to unphysical complex singularities in mesonic propagators. This 
is the case of the proper-time regulator. These complex singularities cause complications. 
For instance, the analytic continuation from the Euclidean to the Minkowski space has 
to be done with care, since the Wick rotation picks up contributions from these poles 
(and also, as we will show, from the infinite semi-circle). Other regulators may be free of 
complex poles {e.g. the Pauli-Villars case) but they have other problems. 

The presence of complex poles is relevant to calculations of hadronic properties in 
effective low-energy models (baryon calculations are reviewed e.g. in Refs. |]I]-§[, and 
meson ones in Refs. [0-|^). Properties of baryons found at the soliton, or mean- field level, 
{e.g. baryon mass, coupling constants, form factors, etc.) are leading-Ai'c results. These 
quantities acquire meson-loop corrections, which are l/Nc suppressed. Since = 3, 
some of these formally suppressed effects may be numerically large. When Wick rotations 
are performed to calculate the meson-loop contributions, the extra complex singularities 
contribute and may not be neglected. 

When discussing the regularization of effective theories, such as the Nambu-Jona- 
Lasinio model, it is useful to remember that regularization is very different from renor- 
malization, in which a cut-off is introduced and later taken to infinity. The problems 
discussed in this paper disappear in the limit of infinite cut-off. In the Nambu-Jona- 
Lasinio model, in order to fit data, we are obliged to work with a finite cut-off which 
rarely exceeds 1 GeV. QCD does not provide us with a prescription for regularizing the 
model. Various regularizations have been used. The proper-time regularization of the 
quark loop is most often used in soliton calculations because it provides us with an action 
which is a local functional of the fields. This simplifies soliton calculations as well as the 
definition of conserved currents, which acquire familiar forms. Proper time regularization 
belongs to the class of regularizations in which only the real part of the Euclidean action 
is regularized. The reason for not regularizing the imaginary part is two-fold. First, the 
proper-time regularization can only regularize a positive definite operator, such as DD"^, 
and it cannot be applied to an operator such as D / D\ the log of which contributes to the 
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imaginary part of the Euclidean action[^. Second, 't Hooft's anomaly matching condition, 
as well as phenomenology, instruct us that the anomalous processes — those which arise 
from the imaginary part of the action — should not be regularized. 

However, other regularizations exist which yield the correct anomalies: four-momentum 
regularizations, for example [|n|-[T^ . These can be expressed either in terms of momentum- 
dependent constituent quark masses (or, equivalently, non-local fields) or in terms of 
cut-off functions which multiply the quark propagators. The non-locality of the resulting 
model yields extra terms contributing to conserved currents. Four-momentum regular- 
izations can and should be applied to both the real and imaginary parts of the action. 
They have the attractive feature that anomalous processes turn out to be unregularized. 
In addition they regularize the theory: both quark and mesons loops are regularized with 
a single regulator. However, the non-locality induced by the four-momentum cut-off reg- 
ularizations make the soliton calculations an order of magnitude longer. That is why the 
proper-time regularization is more often used and why consider it in the present paper. 

A further problem with meson loops is the need for introducing a new cut-off. Indeed, 
the proper-time regularization will regularize the quark loop but not the next-to-leading 
order meson loop corrections. The regularization of the meson loop introduces an extra 
parameter in the model to be fixed by some extra physical requirement. 



2 The regularized action 

We work with the Lagrangian of the SU{2) x SU (2)-symmetric Nambu-Jona-Lasinio 
model with scalar and pseudoscalar interactions in the chiral limit of vanishing current 
quark mass. 



2a2 

The partition function of the system is given by the path integral 



Z = Tt e-^^ = ^j'^^ e-^W , (2) 



2 The imaginary part of the Euclidean action is odd in the e°^^^^ tensor, hence it describes 
anomalous processes. Conversely, the real part, which does not contain the e tensor, describes 



non-anomalous processes. See e.g. [1C|. 
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where is a four-component chiral field, = (0O; 0i; 02, ^s), introduced to partially 
bosonize the model and A/" is a normalization factor. The Euclidean action 5'(0) 



consists of a quark part and a mesonic part, 

2 

^(0) = -Tr logD + ^I^V^ ^ 
with the Dirac operator and the Dirac Hamiltonian defined as 

D = dr + h, h = -la-V + J2Ta^a . (4) 

a 

The quark-meson coupling matrices are Tq = jq and Fj = 27o75Tj, i = 1,2, 3. The trace in 
Eq. d^) runs over the quark field variables. The normalization factor in Eq. (|^) is equal 
to7V = /P0e-^'^'. 

Most of the discussion of this paper concerns the so-called proper-time regularization 
of the action (j^). Indeed it is one of the most widely used regularizations, especially 
in connection with soliton calculations p5|-p!7| . The proper-time regulated action of the 
Nambu-Jona-Lasinio model is 

^a(0) = ^Tr / ^ e-^^'^ + \j dW - ^Tr log(D/Dt) . (5) 

1 

The first term is the real part of the action, which is regularized, whereas the last term is 
the imaginary part of the action, which is finite and is not regularized. In our S't/(2)-case 
with no external sources the contribution of the imaginary part vanishes and we drop it. 



3 Evaluation of the partition function 

We evaluate the partition function (H) using the saddle-point method. The stationary 
field configuration is determined by the equation: 



5S(0) 



. (6) 



We consider the case of spontaneously broken chiral symmetry, where all the fields 0^ 
vanish except one which acquires the non-zero value: 0o = M. We refer to M as the 
constituent quark mass. Performing the Gaussian integration over fluctuations of the field 
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we get: 



logZ = 5(0) + ^Tr \og{K-^/a^) 



where ^ is the inverse meson propagator, defined as 

525(0) 



(7) 



{x,a\K ^ \y,b) 



5(t)a{,x)5(t)i,{y) 



(8) 



The first term in Eq. (|^) is the quark-loop contribution and it is proportional to N^. It 
is the Hartree term. The second term in Eq. (|^ is the meson-loop contribution, and it is 
proportional to A*"^. It contains both the Fock term and the RPA correlation energy. The 
constant under the log results from the norm M . Up to now all of the calculations in 
the soliton sector have been limited to the quark-loop level. In the vacuum sector a few 



attempts have been made to include the meson-loop effects 18-22 



4 Complex poles of the regularized meson propagators 

In this section we analyze the meson propagators K in the vacuum sector when the 
action is regularized by the proper-time method. We show that the propagator has extra 
complex poles, in addition to the expected singularities corresponding to the on-shell 
meson pole and the qq production threshold. This feature has been overlooked in previous 
investigations. 

Since the chiral symmetry of the vacuum is spontaneously broken, all the fields 0a 
vanish except one which acquires the non-zero value: 0o = M. We refer to M as the 
constituent quark mass. From (§) we deduce the following inverse pion and a-meson 
propagators, which are functions of the Euclidean momentum Q: 

K-,\Q^) = ANX,{Q^ + 5ao ^M')f{Q') ■ (9) 

The function /(Q^) is equal to: 

1/A2 ^ ^ 

with Ei{z) = dt/te~^^. Equation (|) implies that the pion propagator (a = 1, 2, 3) has 
a pole at = 0, and that the cr-meson propagator (a = 0) has a pole at = —AM^. 
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These poles are interpreted as on-shell physical particles. On the negative real axis in 
the plane we find, in addition, a quark-antiquark production threshold starting at 

= — 4M^, which is an artifact of the lack of confinement. These are well known 
features of the meson propagators in the Nambu-Jona-Lasinio model. 

The meson propagators have, however, additional singularities in the complex 
plane. They are generated by the exponential regulator which appears in the expression 
( |T0|) of the function fiQ"^). The regulator produces extra zeros of the function / (Q^) in 
the complex plane. This is displayed in Fig. |l|, which shows the analytic structure 
of the pion propagator in the complex plane. We see that, in addition to the 
physical pole and the qq cut, there are infinitely many poles which are located close to the 
imaginary axis. The figure is drawn for the special choice of parameters A = M (which 
corresponds to some choice of the parameter a), but the behavior is qualitatively the same 
for any other choice of parameters. The location of first few poles is: Q^/A^ = 2.2 ± 18.0 z, 
1.4 ± 41. 7z, 0.3 ± 66. 6i, etc. (A = M). Asymptotically, for large | |, the real parts 
of at these poles tend to a constant, and imaginary parts are separated by the value 
87rA2. 

The result described above follows from a simple numerical calculation which can be 
carried out e.g. with Mathematica. We just treat as a complex variable, and evaluate 
the integral over the Ei function in Eq. (|TUp numerically. The function Ei{z) has a cut 
from to — oo, and is analytic elsewhere. The integral in Eq. ( |T0| ) is well behaved, and can 
be evaluated to any desired accuracy. The appearance of complex zeros of the equation 
/(Q^) = is not obvious at first glance, and their location cannot be obtained analytically 
due to the transcendental nature of the function. 

A simple example can show why complex zeros appear in functions of this type. 
Eq. (|l^) involves averaging over exponential functions. Consider a much simpler case of the 
function g{z) = Jq dae"^ = z~^(e^ — 1). This function has zeros at locations Re(z) = 0, 
Im(z) = 27rk (k = 0,1,...), which is very similar to the case shown in Fig. |l[ The function 
of Eq. (0) is clearly more complicated, but the basic structure remains. In Sec. ^ we will 
encounter another example of this behavior. 

The structure such as shown in Fig. |T] is quite peculiar. At first, it might seem that 
causality and unitarity are violated, since singularities appear off the real negative axis in 
Q^. However, as demonstrated by Lee and Wick |2^, and Cutkosky et al. [^, there exist 
models with complex singularities which do satisfy causality and unitarity of the 5'-matrix. 
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This is achieved by restricting the in and out scattering states to be composed of stable 
modes only. Then unitarity of the ^-matrix can be maintained. For a two-point function 
at low values of the s variable, the method of Cutkosky et al. is equivalent to calculate 
Feynman diagrams of the perturbation theory along Wick-rotated contours. At higher 
values of s pinching of singularities may occur, and the prescription is more involved. 
Thus, unitarity and (macro) causality pH] may be satisfied, despite the occurrence of 



complex singularities. An independent study would be required to ascertain whether the 
prescription may be applied in our case. 

Recently Langfeld and Rho |^ analyzed a confining Nambu-Jona-Lasinio-like model 
in which the pion propagator has, in addition to the physical pole, complex singularities 
located similarly to our case. In fact, complex singularities appear frequently in model 
descriptions of confinement (for a review see p7[|). 

A related problem is the asymptotic behavior of the meson propagator for large values 
of I Q P in the complex plane. In derivations of dispersion relations, or in expressions for 
the one-meson-loop energy (see Sec. ^), one makes use of the Cauchy theorem. The unreg- 
ularized meson propagators behave asymptotically as l/Q^, hence the integral over the 
infinite semi-circle vanishes (after a sufficient number of subtractions), and one obtains 
the usual dispersion relations. With the finite proper-time cut-off, the asymptotic form of 
the meson propagator has a form involving an exponential function, ~ exp(— ^)/i((5^), 
where /i is a slowly varying function. In applications of the Cauchy theorem, the contribu- 
tion from the infinite semi-circle does not vanish. This problem will be discussed in more 
detail in Sec. ||. 

We summarize the problems encountered in the Nambu-Jona-Lasinio model with a 
proper-time regulator: 

(1) The meson propagators have (infinitely many) extra poles which lie close to the imag- 
inary axis in the complex plane. After performing the 3-momentum integration, 
these poles unwind in cuts in the energy variable. 

(2) When performing a Wick rotation, the complex poles contribute. In addition, we get 
contributions from the infinite semi-circle. As a result, the usual dispersion relations 
for meson correlators and other Green's functions acquire extra contributions. 

(3) Note, however, that the pathological poles are quite far away from the origin in the 
complex plane. For the case of Fig. |I| (where M = A) the closest poles are about 
18 units of A^ away, while the qq production threshold is 4 units of A^ away. Thus 
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in some numerical calculations the presence of these poles may not be immediately 
noticeable. 

(4) In calculations of Green's functions at large values of incoming momenta squared, a 



prescription like the one in Ref. has to be used in order to satisfy basic require- 
ments of field theory. 

All above problems are due to the finiteness of the cut-off. In the limit A ^ oo the poles 
at complex locations are moved to infinity, and the usual analytic structure is recovered. 



5 Spectral representation of meson propagators in background soliton fields 

In this and the following sections, we limit the number of quark states contributing 



to the spectral sum ([T^ ) so as to maintain the propagator K finite even in the limit of 
infinite cut-off. This allows us to obtain a better understanding of the extra poles which 
are produced by the proper-time regulator. 

We extend our results to the general case with arbitrary background fields, e.g. for 
solitons. We make use of the spectral representation, given by the eigenstates of the Dirac 
Hamiltonian with the meson fields in a saddle-point configuration: 

^ I j) = (^-ia + 12 Ta^a^ I j) = I j) , (X I j) = qj{x) . (11) 

A straightforward calculation, outlined in App. 0, yields the following expression for the 
inverse meson propagator 

(f, a I K-^{y) \y,h) = 5{x - y)5aba^ + [iW)^ aQkix)) (ql{y)T bqj{y)) x 



jk 



ds / _„,2 „,2 / , ,a o\ f du 



1/A2 ^ \ -1 



(12) 



Note that the sum over the indices j and k runs over all the quark states and that it 
includes a sum over color. In the infinite cut-off case Eq. (^) reduces to the standard 
expression for the RPA propagator [p8H 30|l: 
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(f , a I K^^^{u) \y,b) =5{x- y)5aba? - 

a 

^\ ep-th + iv Cp - th - iv 

where the labels p and h denote empty particle (e^ > 0) and occupied hole (e/^ < 0) 
states. (A finite quark chemical potential /ig can be introduced so as to include possible 
valence states in the set of occupied hole states.) In expression (|I^) the sum is restricted 
to particle-hole excitations. This follows "automatically" from Eq. ([T^). As A ^ oo, the 
contributions to the spectral sum from pairs of states with the same sign of the energy 
vanish identically. This the expected Pauli-blocking effect. 




6 Energy of the soliton at the one-meson-loop level 

For the case where the background fields are stationary, the expression for the energy 
of the system is 

E = -^\ogZ = EQ + Ei- vac , (14) 

where Eq is the contribution obtained in the saddle-point approximation containing the 
quark loop, and Ei is the contribution from the one-meson loop. The vacuum energy (vac) 
is subtracted. The leading saddle-point contribution, which is of order Ac is: 

^0 = i^E/ + y/ t'x^iJ^f ■ (15) 

All soliton calculations up to now have been limited to this leading-order contribution. 
In order to make contact with familiar many-body theory, we define the quark-quark 
interaction: 




Making use of Eq. (|]) we can cast the quark-loop contribution in the form: 

^ I ^ [ ds _g,2 Iv^T^ 1 f as _^,2 r ds _^,^2 

= E y 5575^ ' - 2 ? ''"•"S / in''' ' i 7175*' ' ■ '''' 

^ J 1 J« 1 1 

The second term is recognized as the Hartree energy shown diagrammatically in Fig. |^. 
In the limit A — >^ cxd, the expression (|l^ reduces to the familiar expression of the Hartree 
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energy: 



lim Eo = y2^h - ]:y] Vhh,h'h' 



The one-meson-loop contribution is 



oo 



E. = \j ^Tr log [K-\u)/a') =\j£^og Det [K-\u)/a^) , (19) 

— oo — oo 

with given by Eq. (p!2D . In the infinite cut-off limit we may use the Cauchy theorem 
to evaluate the integral over u in Eq. (|19D . As shown in App. the integral of a log of 
a (well behaved) function / is proportional to the sum over all zeros of / in the upper 
complex plane, minus sum over poles of / in the upper complex plane. In our case / 
equals to Det{K][^^/a'^). The zeros correspond to the eigen frequencies u which satisfy 
the equation 

BetKxi^{u) = , (20) 

which is the usual RPA equation (see also App. 0). The poles of DetK]^\^^ are located 
at the particle-hole excitation energies, as seen from Eq. (0). As a result, we obtain 

^im E, = -j £\og Det {k-^1^{v) / a^) = - E - 3 ' ' (21) 

_oo ijJi>0 ph 

where 

cui = Im{u^) (22) 
In App. y we show that Eq. ( PT| ) may be rewritten as 

lim El = --yuilYil^ -^ + -yVph hp- (23) 

We recognize here the usual RPA expression, with a zero-mode contribution |28-0O|]- The 
last term is the Fock term, represented diagrammatically in Fig. ^(a). The ring diagrams 
of Fig. §(b-c) represent the RPA correlation energy. Note that the meson-loop energy is 
not just the sum over the RPA frequencies of the form ^J^i^i^ as sometimes claimed 



3T[| . The form of Eq. (p3|) comes from the qq substructure of the meson propagator. 



Indeed, in purely mesonic models such as the Skyrme model, Ei = at the iV° 

level. 
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As mentioned before, Ei is suppressed by one power of Nc compared to Eq. However, 
many-body physics provides examples of 1/A^c-suppressed results which are important. 
One effect is the so-called zero-mode correction, which accounts for spurious contributions 
of the center-of-mass motion. The mean field approximation used in soliton calculations 
does not separate the center-of-mass coordinates from intrinsic coordinates. This notorious 
problem of projecting out the center-of-mass motion of solitons has attracted a lot of 



attention |l31|-p8| and numerical estimates give a negative contribution of the order of 300 
MeV. Some estimates use an expression derived for systems with a finite number of 
particles in the RPA approximation, as in Eq. (p3D 



AM = — (24) 
2M ^ ' 

where P is the momentum operator, and M the soliton mass. In general, such a zero- 
mode subtraction appears for every continuous symmetry of the Lagrangian which is 
broken by the solitonic solution. In addition to the translational invariance, hedgehog 
solitons appearing in chiral models also break the rotational and isospin symmetry, leaving 
only the sum of spin and isospin preserved. In analogy to Eq. ( p4|) this leads to another 
subtraction AM = — (J^)/(20), where J is the angular momentum operator, and B is 
the moment of inertia For typical model parameters, the rough estimates of the 
zero-mode corrections are large, up to 30% of the soliton mass, as expected. 

In the derivation of Eq. (^3]) the Wick rotation has been used. As will be shown in 
next sections, extra contributions arise when the the proper-time regulator is present, and 
expressions as Eq. (|23| ) do not hold. 



7 Finite proper-time cut-ofF 

We decompose the integrals over s as follows: 



oo oo 1/A 

J ds = J ds - J ds , (25) 

1/A2 

and use 



' ' 
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where the error function erf (2;) and the incomplete Euler 7 functions 7(n, z) and 7*(n, z) 
are defined e.g. in Ref. We decompose the total contribution to into the infinite 
cut-off part, and a remainder, which describes the finite-cut-off effects: 

K-\u) = Kxl^iu) + AK-\u) , (27) 

where 

(x,a I AK-\u) I y,b) = 1 ^ (g](f)r,g,(f)) (gi(y)r,g,(y)) x 
erf(M) ^ erf(^) \ ^ 

4A^ y^"^ U' ^ 2A^^ 

-1 ^ 

The function 7* (a, 2;) is an entire function of the variable z {i.e. it is single-valued and has 
no singularities in the finite complex plane - singularities occur only at infinity). Therefore 
{x, a I AK^^lu) I y, h) is also an entire function in the variable v. This is because v enters 
in the argument of 7* in a non-singular combination with the integration variable u. In 
other words, AK~^{v) has no singularities in the finite complex plane. As a result, the 
singularities of the function K~^{v) are generated by its A 00 part in Eq. (pT]). This 
means they occur at particle-hole excitations, v = ±i(ep — e^). 

The zeros of DetK^^^u) will of course be at different locations than zeros of the 
A ^ 00 part. As we change A from infinity to a finite value, the zeros of Deti^'"^ move 
away from those of DetK][^^. In addition, new unphysical zeros emerge at complex values 
of u. The situation is analogous to the case of the vacuum meson propagators, discussed 
in Sec. |[ 





8 Two-level model 

To illustrate how the new zeros of Det-ft' are generated, we consider a two-level model, 
in which just one particle state p and one hole state h contribute to the spectral decom- 
position (0) of the propagator K. In this case only one mesonic fiuctuation occurs. Fig. || 
shows the analytic structure of the corresponding inverse meson propagator K^^ (z/) in 
the complex plane of the energy variable u. We note the poles (filled dots) located at 
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±z(ep — eh). The physical zeros (empty dots) on the imaginary axis correspond to the 
RPA vibration energies. We also see infinitely many unphysical zeros (empty dots) close 
to the axes at 45° and 135°. If it were drawn in the z/^ complex plane, Fig. ^ would look 
very similar to Fig. |I]. The displayed structure results from a numerical calculation, where 
we evaluate expression (|T^ for the two-level case. 

We can carry out an even simpler analysis. Notice that for large values of | z/ |, the 
function in the two-level model can be approximated by: 

K-\u^)/a^ c:^l + b{u^)e^ , (29) 

where &(i^^) is a slowly varying function. Assuming &(z^^) to be a real positive constant, 
and breaking up z/^ into real and imaginary parts, z/^ = + il, the zeros of the function 
occur at the locations 

/ = 4A2(2fc + l)7r , R = -4AHogb , (30) 

where k is an integer. The difference between the imaginary parts of successive zeros is 
SttA^. If we integrate by parts, as in Eq. ( |B.1| ), then the integral to be considered is 

J= j du\og{K-\u)/a') = - J Vv " / '^''^^''^ ■ ^^^^ 

— oo — oo — oo 

Denoting u =\ u \ e^'^ , we find that for large | u \ 
~ for cos(20) > , 



i(^) ~ ^^''^ for cos(20) < . (32) 

This allows us to check if we can close the contour in order to evaluate J via the Cauchy 
theorem. As a consequence of Eq. (P^D, the integral of j(z/) along the contours Ci and C3 
in Fig. ^ vanishes (this is the case cos(20) > 0), but the integral along the contour C2 
(case cos(2</)) < 0) does not vanish, and is equal to V2 \ v |^ /(6A^). Hence the integral 
of j over the contour C2 goes as | v |^. 

The result is that J in Eq. (|3TD is an infinite sum of the pole contributions, plus an 
infinite contribution from the contour C2. Since we know that J is finite, cancelations 
between the infinities from the pole and the parts occur, and we are left with a finite 
(but non-zero) result. Certainly, the same cancelation of infinities takes place in the general 
case of Eq. (pj]), when we evaluate it via the Cauchy theorem. 



|2 
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9 Other regulators 



The analytic structure displayed in Figs. |l| and || is specific to the proper-time regula- 
tor. Other popular regulators, such as the sharp Euclidean four-momentum cut-off or the 
Pauli-Villars regulator do not lead to complex singularities. The Pauli-Villars regulariza- 
tion consist of introducing additional families of quarks (which have Bose statistics). In 
the semi-bosonized form the action has the form 

C = -Tr log(-27 • d + -foTa(pa) + CpTr log(-i7 ■ d + gp-foTa(pa) + ^ / d'^^4>'^ >(33) 

where the the label p labels the Pauli-Villars quark families, is the number of these 
families, and the constants Cp and are chosen appropriately in order to cancel infinities 
4^. It is clear from Eq. ( pSj) that the corresponding meson propagators will have singu- 



larities along the real negative axis only. Apart from poles, there is the cut going from 
the qq production threshold at = — 4M^ to = — oo, and additional cuts associated 
with the production of the Pauli-Villars quarks, going from = —Ag^M^ to = — oo. 
Therefore the Nambu-Jona-Lasinio model with the Pauli-Villars regulator obeys causality. 
Moreover, the usual analyticity structure allows for the straightforward use of the Cauchy 
theorem, rotation of contours, etc. In particular, this leads to the formal expression for the 
RPA correlation energy, such as in Eq. (p3D . The problem with the Pauli-Villars regulator 
is that, as remarked in Ref. it violates unitarity. This is because for sufficiently large 
negative values of {—Q^ > 4:{M'^ + ^f^M^)) the discontinuity along the cut in meson 
correlators is negative. 



10 Conclusion 

We have studied the meson propagators in the Nambu-Jona-Lasinio model in the 
vacuum and solitonic backgrounds. We have found that the proper-time regularization 
introduces extra complex poles to the meson propagators, which may cause serious com- 
plications. The extra poles, as well as the non-vanishing contribution of the infinite semi- 
circle, forbid the usual deformation of the energy integration contour which allows one 
to express the partition function in terms of the physical modes of excitation of the sys- 
tem. The unphysical complex singularities, and the infinite semi-circle piece, have to be 
accounted for. As a result, the formal connection to well-known text-book expressions is 
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lost. We no longer have the usual dispersion relations for Green's functions, the RPA- 
equations for the collective excitations of the system can not be derived in the usual 
symplectic form. In general, we cannot use the particle-hole description. 

In view of this the naive use of the Cauchy theorem in order to calculate the RPA- 
correlation energy for the Nambu-Jona-Lasinio soliton as e.g. in Ref. is incorrect. 

The correct way to calculate the energy and other observables on the one-meson-loop 
level is to perform explicitly the integrals over Euclidean four-momenta in the meson 
propagators. 

Our remarks concerning meson correlators also apply to diquark correlators. Diquark 
models have been formulated in the spirit of the Nambu-Jona-Lasinio model [^4|-|49[|. 
Our warning is that similar care has to be taken while performing integrals over diquark 
momenta, e.g. when continuing diquark functions from the Euclidean to the Minkowski 
space, or solving the Fadeev equations in quark-diquark systems. 

The authors acknowledge the support of the Polish State Committee of Scientific 
Research, grant 2 P03B 188 09 (WB), the Maria Sklodowska-Curie grant PAA/NSF-95- 
158, the Bulgarian National Science Foundation, contract $-406 (EN), and the Alexander 
von Humboldt Stiftung (GR, WB), which made this collaboration possible. 
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A Meson propagators in the spectral representation 



We expand the first term in Eq. to second order in mesonic fluctuations Scf) around 
the reference stationary state (vacuum, soliton), which has the meson field configuration 
(j). We decompose the D'^ D operator as follows: D'^D = Hq + Vi + V2, where 

Ho = -dl + h^ , 

a 

F2 = Era50aWb . (A.l) 
ab 

Using the expansion 

1 } 



we pick-up the second-order piece of the proper-time effective action expanded in the 
mesonic fluctuations 50: 

SF{2) ^ J Ly^e-^^o + dpVre-'^^-^^'''^V,e-'^A . (A.3) 

1/A2 \ / 

Now let us use the spectral representation, with states denoted by | uj,i), etc., where oj 
is the frequency variable (conjugate to r) and j is as in Eq. |TT|. Note that Hq is diagonal 
in this representation: 

Ho\uj,j) = {uj^ + t])\u,j) . (A.4) 

We obtain 

(f , a I K-\v) I h) = Six - y)5a,a' + (S^ + T^^) , (A.5) 

jk 



where 



Sil = -\ [i{x)V,q^{x)) (gi(y)r,g, (y)) / ^ J rf.e-(-^+^?) + (a ^ 6, f ^ f) , 

1/A2 

^ (g](f)r,gfc(f)) (gi(y)r,g,(y)) | ^ ((e, + e^f + z/^) 



rpik 



oo 1 



X J dssj dpe"^'-"'^^^-^' +^fe^e-*^^" + {a^b,x^y) . (A.6) 

1/A2 
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Carrying out the a; integration, changing the integration variable /3 to u = 2(3 — 1, and 
noticing that the replacement {a ^ b,x ^ y) is equivalent to (j ^ k, u ^ —u), we get 
Eq. (0). 



B Integral of the logarithm of a function. 

We recall a useful form of an integral of the logarithm of a function /. First, integrate 
by parts: 

OO OO 



00 
—00 



Assume that the function / has poles at z/ = z/f°''' and zeros at u = Uj. Then the function 
^ '^^'^fil)" ^ poles at u = uf"^^ with residues —1/?°'*^, and another set of poles 

at u = Uj with residues If the surface term vanishes, and if the integral along the 
infinite upper (or lower) semi-circle vanishes, then the Cauchy theorem gives the following 
result: 



'J ~ ' 

U>0 i>0 



where j > indicates that Xm(z/°) > 0, and i > indicates that lm{h'f°^'^) > 0. 



C RPA equations. 

It is convenient to use an orthonormal basis in the mesonic space, defined by a set of 
functions fn{x). Each /„ is a four-component vector in the a, tti, 7T2, and tts space, i.e. 
/n = (/°,/^,/n',/n')- We define 

K-'niiy) = T.Jd'xJ a I K-\u) \ y, b)f:{y) . (C.l) 

ab 

Condition (pO]) is equivalent to solving the eigenvalue problem (everywhere below repeated 
indices are summed over) 

hm K-'Mf^ = . (C.2) 

A^oo 
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Defining the quark-meson overlap coefficients as 



(C.3) 



we can show that the problem (|(J.2D is equivalent to the following set of equations: 



(ep-e, + zz/)X^"^ + (Cf)*/n = 



(C.4) 



which can be seen immediately when one eliminates the variables X^'^ and Y^'p from the 
above equation — as a result, the original Eq. ( |C.2| ) is obtained. One can also eliminate the 
variables /„ from Eqs. (^), using the first of Eqs. (|g: /„ = -l/a^ (Cf Xp'^ + C^pY^p). 
The second and third equations become: 

(ep -eh- iu)Y''P - \ ((C^^')*C^''^'XP'^' + (C^^')*C^'pV^'^'') = . (C.5) 

Since the basis {/«} is complete (this is an arbitrary orthonormal basis, e.g. of plane 
waves), i.e. J2n fni^)fniy) = ^''^^{x-y), we have 



ij,lk ) 



(C.6) 



where V has been defined in Eq. ([I6|). Using this notation we may rewrite Eq. ( p.5|) as 



(Cp — ^h)Sph,p'h' + Vph,h'i 



Vhp,h'p' 



Vph,p'h' 



\ ( \ 

J^p'h 



(cp — ^h)Sph,h'p' + Vhp,p'h' j yY'^'p J 

( \ ( \ 

Sph,p'h' 



—Sph,h'p' 



yh'p' 

V " J 



(C.7) 
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This is the symplectic RPA eigenvalue problem, familiar from nuclear physics, which has 
the form 



A B 



\ 



B* A* 



J 



( \ 

X 



I 



—IV 



\ 







1 



X 



(C.8) 



We refer the reader to Refs. [Pq-pO| for a detailed discussion of this symplectic eigenvalue 
problem. We now want to rewrite the quantity J2phi^p~^h) differently. From Eq. (|C.7[ - P.8|) 
we see that 



^(cp - e,,) =TtA-J2 Kh,hp 

ph ph 



(C.9) 



The trace in this expression may be evaluated using the eigenstates of Eq. (|C8|) , and we 

get 



«>0 



(C.IO) 



We stress that the above result is algebraic in origin, and holds for any symplectic problem. 
The zero-mode piece ^m' is generated "automatically" . In hedgehog models there is also 
an angular-momentum zero-mode piece. 

Collecting expressions ( pT| , |C.9| , p.lO|) and using the fact that | Xj ^ — | p= 1 
2|-g, we get finally 



1 (P^) 1 

lim Ei = — y^cOilYil^ + - y^Vph hp 

i>0 ph 



(C.ll) 
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Fig. 1. The analytic structure of the pion propagator in the complex Euclidean momentum plane, 
Q^. It has the physical pole at = 0, the qq production cut from = — 4M^ to — oo, and 
unphysical poles located close to the imaginary axis. The parameters are M = A. Note different 
scales on the real and imaginary axis. The corresponding plot for the cr-meson propagator differs 
only by the location of the physical pole, which is moved from = to = — 4M^. 
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Fig. 2. The Hartree diagram. The dotted hne is the interaction, which is local in our case, and 
carries a factor oil/Nc- The diagram is of order Nc- 




(C) 

Fig. 3. The Fock digram (a), and higher order ring diagrams (b-c). All these diagrams are of 
order N^. 
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Fig. 4. The analytic structure of the function K~^{i>). It has poles on the imaginary axis (filled 
dots), corresponding to particle-hole excitations, and zeros on the imaginary axis (empty dots), 
corresponding to the RPA vibrations. Additionally, it has infinitely many unphysical zeros in 
the complex plane off the imaginary axis (empty dots). These zeros become poles of the meson 
propagator K. The contour integral along C3 diverges as the radius of the semi-circle goes to 
infinity. The sum of pole contributions and the Ca-contribution is finite. 
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